TeV Neutrinos and GeV Photons from Shock Breakout in Supernovae 
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We show that as a Type II supernova shock breaks out of its progenitor star, it becomes coUisionless 
and may accelerate protons to energies > 10 TeV. Inelastic nuclear collisions of these protons produce 
a ~ 1 hr long flash of TeV neutrinos and 10 GeV photons, about 10 hr after the thermal (10 MeV) 
neutrino burst from the cooling neutron star. A Galactic supernova in a red supergiant star would 
produce a photon and neutrino flux of ~ 10"* erg cm~^ s~^. A km^ neutrino detector will detect 
~ 100 muons, thus allowing to constrain both supernova models and neutrino properties. 

PACS numbers; 97.60.Bw, 14.60. Pq, 98.70.Rz, 95.85.Ry 



Type II supernovae are triggered by core collapse in a 
massive star, generating a strong shock wave which prop- 
agates through the progenitor star and ejects its envelope 
(see for a review). Interest in the breakout of the 
shock out of the envelope of the star, was prompted by 
Colgate's suggestion that it may be accompanied by ener- 
getic 7-ray radiation due to Brenisstrahlung and inverse- 
Compton emission of thermal electrons ||^. Numerical 
simulations concluded that the shock breakout should be 
accompanied by a burst of hard UV and X-ray radiation 



In the optically-thick interior of the progenitor, the 
shock is dominated by radiation. The dissipation of ki- 
netic energy at the shock front is mediated by Comp- 
ton scattering rather than by collisional viscosity, and 
the thermal energy density behind the shock is mostly 
in the form of radiation. Numerical simulations suggest 
that during breakout, the radiation-dominated shock is 
transformed into a viscous shock ||,^ . When the optical 
depth of the gas lying ahead of the shock is no longer very 
large, radiation decouples from the gas and the shock is 
transformed into a collisional, viscosity-mediated shock. 

In this Letter, we show that as the shock becomes 
viscous, electro-magnetic instabilities at the shock front 
grow at a rate which is many orders of magnitude higher 
than the collision rate, leading to scattering of ions and to 
magnetic field amplification on a time scale much shorter 
than the viscous time scale. Under similar conditions, as- 
trophysical shocks are known to become coUisionless (i.e., 
to be mediated by collective plasma instabilities rather 
than by collisional viscosity) and to accelerate charged 
particles to relativistic energies (see ^ for a review) . We 
derive the expected flux and spectrum of high energy 
photons and neutrinos which are produced through in- 
elastic nuclear collisions of the accelerated protons. 

As we show below, the resolution of available numeri- 
cal simulations is not sufficient to correctly describe the 
evolution of plasma parameters across the viscous shock 
(see also [^). We therefore first derive the expected col- 
lisional shock structure from basic principles, and only 



then discuss the stability of the shock. While our analy- 
sis provides a description of the shock structure, it does 
not allow to determine the exact radius where the vis- 
cous shock forms. Numerical calculations suggest that 
the transition to a viscosity-dominated shock occurs at 
an optical depth (measured from the stellar surface) be- 
tween r ~ 1 1^ and r ~ 10 0. We therefore parame- 
terize our results in terms of r. More detailed numerical 
simulations are necessary to allow an accurate determi- 
nation of this parameter. 

We focus our discussion on stars which are believed to 
be the most common progenitors of Type II supernovae, 
namely red-supergiants (RSG) [Q. Our model predicts 
emission of high energy particles also for blue-supergiants 
(BSGs), such as the unusual progenitor of SN1987A. 
However, the predicted neutrino and 7-ray fluxes from 
the surface layer of BSGs are predicted to be two orders 
of magnitude smaller than from RSGs. Additional emis- 
sion may result from the interaction of the supernova 
shock with the progenitor wind on long time scales 
although this emission depends sensitively on the wind 
speed and mass loss rate and might also be strongly mod- 
ulated by past variability of the wind . 

Physical parameters at shock breakout. When the 
radiation-dominated shock reaches an optical depth lower 
than unity, radiation can no longer couple to the gas and 
so the shock becomes viscous. This occurs very close to 
the surface radius of the progenitor star, Ri,. For the 
high temperatures of interest, the opacity is dominated 
by Thomson scattering, and the mass of a shell of optical 
depth T located ahead of the shock is given by 
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where i?^.i3.5 = (i?*/10 cm) = 1 reflects the charac- 
teristic radius of a RSG progenitor. As mentioned above, 
we ignore rarer, more compact progenitors, such as blue 
supergiants for which i?^,i3.5 = 0.1, since the shell mass 
at a fixed r, and hence the energy output in high-energy 
neutrinos and photons, are smaller by a factor 10^ in 
them compared to RSGs. Using the analytic approxi- 
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mations of Matzner & Mckee Q to describe the outer 
envelope structure, we find that the mass density of gas 
at the RSG surface is 



p = 0.54 X lO^^V/^ 



g cm 
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and the shell thickness is A = 1.1 x IO^^t^^^ cm. 

The shock velocity at breakout is Vg ~ lO^cm s^^ (e.i 
PI ) , corresponding to an ion temperature of 



T, 
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1.0 X lO^u^Q eV, 



(3) 



where fj, = mp/2 is the mean atomic mass per particle for 
a plasma in which the thermal energy is shared equally 
between protons and electrons. If the shock is radiation- 
dominated, the radiation temperature is approximately 
given by the relation aT^ « pw^, i.e. 
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eV, 
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(p/10"i° g cm-3) 



where p-io 

Viscous shock structure. The ion-ion, ion-electron, 
and electron-electron collision rates are given by 
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in3 fsr, 5^^^S \ if Te < TOeT,/mp; 

otherwise, 



lOV-ioT, 



-3/2 



e,4 



(7) 



Throughout this Letter.^ we use the notation T^.x ~ 
(Ta/lO^^ eV) to denote the temperature of component 
a (= i for ions, e for electrons, or r for radiation) of the 
plasma. The collisional energy loss of ions is dominated 
by ion-ion collisions once the electrons are heated to a 
temperature > 0.07T,. 

The electrons cool by Compton scattering and 
Bremsstrahlung emission. Compton scattering domi- 
nates for Te > 10^ eV, and provides a cooling rate of 



l^Comp 



iaT/SmeC)aT^ = ST^^ g s'^. The electron 



temperature is determined by the balance between heat- 
ing and cooling, VieTi w z^compTe, yielding a post-shock 
electron temperature 
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3 X lO\T,,5P-wy/'T-li' eV. 
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At this temperature, the ion-electron collision rate is 



(9) 



Based on Eqs. we conclude that the ion ki- 

netic energy dissipation is dominated by ion-ion colli- 
sions, leading to a shock front width Ag^u w v'l^Va — 
3 X lO^plJo The electrons are heated by the ions to 



a temperature Te+ « 30 keV which is lower by a factor 
of a few than the ion temperature, Ti sa 100 keV. Behind 
the shock, the electrons cool by Compton scattering, but 
their temperature is maintained at Tg « Te-i- over a shell 
of width « v~}vs ~ 10®p_Jq^ cm, due to their coupling 
to the ions. Over a larger length scale, both the electrons 
and ions cool below their post shock temperature. The 
electrons maintain a local thermodynamic equilibrium, 
since they are heated to Tg ^ O.OTTi over a length scale 
<C As.ii and subsequently Vee ^ t'ie- 

Our results agree with the numerical simulations of 
Klein & Chevalier ^ for a RSG progenitor, who con- 
cluded that a collisional shock is formed at t ~ 1 where 
Tr ~ 10 eV. Klein & Chevalier derived a post-shock tem- 
perature of ~ 10** eV, but pointed out that this result 
depends on the artificial viscosity they used, with higher 
temperatures reached for lower artificial viscosity values 
(which yield a sharper shock), and thus concluded that 
the gas should be heated up in reality to T « lO^eV. 
Ensman & Burrows |^ found in their simulations of a 
BSC progenitor that a collisional shock forms at r > 10 
but reaches lower temperatures, ~ 1 keV. This calcula- 
tion, however, assumed thermal equilibrium between the 
ions and the electrons. Unfortunately, the resolution of 
all numerical simulations so far has not been sufficient to 
correctly determine the balance between collisional heat- 
ing and radiative cooling on the relevant length scale of 

10** cm. 

Electro-magnetic stability and the collisionless shock. 
The ion velocity distribution within the collisional shock 
front is highly anisotropic over a length scale A^^^. This 
anisotropy leads to the development of electro-magnetic 
instabilities. In order to analyze the shock front sta- 
bility, we consider a plasma composed of electrons in 
thermal equilibrium at a temperature T^, and protons 
with a bi-Maxwellian distribution having a temperature 
T± = mpVg along one axis (representing the velocity dis- 
persion associated with the ion streams perpendicular to 
the shock front), and a lower temperature in the orthogo- 
nal plane Ti < T± (parallel to the shock front). Since the 
electrons in our case have large thermal speeds, namely 
Te/T± ^ 0.1 3> me/rrip, electrostatic instabilities such 
as the two-stream instability, are suppressed. We there- 
fore focus on electro-magnetic instabilities driven by the 
ion anisotropy. Such instabilities are typically expected 
to grow at a rate of < {vs/c)vpi, where Vpi is the ion 
plasma frequency (e.g. |^), with the maximum rate ob- 
tained only in the limit of infinite T^. 

Solving the dispersion relation for electro-magnetic 
waves in the above plasma, we find purely growing modes 
(i.e. modes having purely imaginary frequencies) with an 
electric field along the direction of the high ion temper- 
ature (i.e. perpendicular to the shock front), as long as 
the wave vector k satisfies the condition ikcjl'KVpi)^ < 
(T±/Ti) — 1. In the limit me/rrip <C Te/Tj_ <^ mp/mg, 
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the growth rates are 



ble or larger than the proton Larmor radius. The Lar- 
mor radius of the thermal protons in an equipartition 



vem 



in(l,fc/F / fc*) , for T^jTi ^ Y, Held is comparable to the wavelength ci^p/ over which 



2^ - Ij fc/r3/2fc,, forTx/T, <r, 



(10) 

where fc* = 2'Ki>pi/c and Y ~ {mpTj_/meTe)^/^ . 

For large anisotropy, T±/Ti ^ Y , the maximal growth 
rate is {vs/c)vpi. If the anisotropy is smaller, the maximal 
growth rate is {vslc)vpi/Y^/'^, which in our case yields 
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The growth rate of electro-magnetic instabilities is there- 
fore many orders of magnitude larger than the collision 
rates va, vip. Correspondingly, the characteristic length 

1 1/2 

scale of the fastest growing modes, cv^^ = 20p_io ^m, 
is many orders of magnitude smaller than the coUi- 
sional shock width. Under such conditions, we expect 
the shock to become coUisionless, i.e. to be mediated 
by plasma instabilities rather than by collisions. The 
Ohmic dissipation rate of the electro- magnetic fields, < 
{nie / rap){i / ci>~^)~'^ i>ee^ is much lower than their growth 
rate for all relevant temperatures and length scales i ^ 

If a coUisionless shock indeed forms, it would accelerate 
charged particles to high energies, producing a power-low 
distribution of particle number, dN /dp oc where p 
is the particle momentum |^]. The acceleration rate of 
either relativistic or non-relativistic protons is given by, 
i^acc ~ i^siE/ Eth)^^ , where E is the kinetic energy of the 
accelerated proton, Eth is the energy of thermal protons, 
and vb is the (non-relativistic) proton gyro-frequency, 
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i^B = 7^ = 1.5 X lO^Bi 
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where Bi = {B/10 G) is the magnetic field normal- 
ized to 10 G, which is the typical value inferred at the 
surface of RSGs |]lO|. Note, however, that the electro- 
magnetic instabilities build-up the magnetic field. In 
analogy with coUisionless shocks in supernova remnants 
| pT[ , 7-ray bursts [Q, or the intergalactic medium [ p^ , 
the magnetic field may grow up to a sizeable fraction of 
the equipartition value, -Bo. p. = lO^'^pj^^QWs.g G, and sat- 
urate near this value due to nonlinear effects Thus, 
the acceleration rate is likely to significantly exceed the 
rate implied by Eq. (p^. Since the acceleration rate at 
thermal energies, vb, is much higher than the collision 
frequency of thermal protons, and since the coUisional 
energy loss rate decreases as E~^^'^ while the accelera- 
tion rate declines only as E~^, protons are expected to 
be accelerated to energies well above thermal. 

Scattering of protons by magnetic field fluctuations is 
most effective if the fluctuation wavelength is compara- 



instabilities grow fastest. As the Larmor radius of highly 
relativistic protons is much larger, a conservative lower 
limit to the maximum energy of the accelerated protons 
is given by 



^max ~ -eBA = ISvs.qBit^/^ TeV, 



(13) 



where we have assumed that instabilities do not lead to 
signiflcant amplification of the field on large scales. This 
is probably an underestimate of the maximum energy, 
as nonlinear effects are expected to amplify the field on 
scales :s> cv~^ . If this occurs, then protons would be 
accelerated up to energies well in excess of 10 TeV. 

Neutrinos and photons. Protons accelerated to rel- 
ativistic energy lose energy by inelastic nuclear colli- 
sions. The average number of nuclear collisions encoun- 
tered by a proton is Nc ~ nuQcA/vg = 5.2t/vs^q, where 
(To = 50 mb. Since a proton loses ^ 20% of its energy in 
each collision, the total energy lost to pion production is 

E^ « ^pmin(l,T)— = 3epmin(T,T2) x lO^^erg, 

TUp 

(14) 

where ^p is the fraction of the post-shock energy density 
which is converted to relativistic protons (note that for 
a spectrum dN/dp cx p^^, the energy of super-thermal 
particles is dominated by relativistic particles). For the 
strong shocks around supernova remnants, the inferred 
values of are of order unity ||l5[| . 

Roughly a third of E^ is converted to muon neutrinos, 
and a similar fraction to high energy photons, through 
pion decay. The neutrino and photon signals should be 
coincident and spread over a time scale of 2i?*/c « 1 
hour. (The emission from the back side of the progenitor 
star is visible for TeV neutrinos but not so for 7-rays, and 
consequently the 7-ray fiash should be shorter by about a 
factor of two than the neutrino pulse.) The predicted flux 
of 7-rays and muon neutrinos from a Type II supernova 
in our Galaxy is therefore 



« F^^ « 10-%min(T,T2)d: 



lolpc erg s ^cm ^ 



(15) 



where diokpc is the supernova distance in units of 10 kpc. 
The high density of ^ 10 eV photons from the radiative 
shock, will lead to a high pair-production optical depth 
for photons with energies > 10 GeV. Thus, we expect a 
photon spectrum dN-^/dE-^ oc E^'^ over the energy range 
lOOMeV < E^< lOGeV. This signal is easily detectable 
with planned 7-ray telescopes (e.g. jlBf). 

The neutrino spectrum is expected to follow the ac- 
celerated proton spectrum, with a number per energy 
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interval of dN/dE oc E~'^ at relativistic proton ener- 
gies. The probability that a muon neutrino will pro- 
duce a high-energy muon in a terrestrial detector is 



1.3 X 10"^^;: 



i/,TcV ' 



with /3 = 2 for i?j/,ToV < 1 and 



13 = 1 for -Bj/^TcV > 1- Thus, the number flux of muon 
induced neutrinos can be related to the neutrino energy 
flux, — J dE^E^{dJv/dEv), where J„ is the neutrino 
number flux, through the relation 



1 + ln(£:,,,nax/lTcV) Pq 
ln(^.,max/lGcV) E^ 



(16) 



where Pq/Eq = 1.3 x lO^^TeV"^ Using Eq. (H) we 
find that 



0.15 



Pq E^/3 
En 47rd2 



130epmin(r,r2)d-otpckm-' (17) 



muons are expected to be detected for a single super- 
nova. The factor 0.15 follows from Eq. (|l^) under the 
assumption i?i/,max ~ 1 TcV. This is probably a conser- 
vative estimate, as proton energies Ep 3> 10 TeV may be 
achieved if the magnetic field is amplified. 

Implications. A type II supernova in our Galaxy, 
which should occur once per 40 years ||l8|, is expected 
to produce a strong neutrino signal of ^ 100 events at 
^ 1 TeV spread over ^ 1 hr, in planned square-kilometer 
neutrino detectors (see, e.g., ||l^ for a review). This is 
the strongest astrophysical signal predicted for the next 
generation of TeV neutrino detectors. Detection of the 
neutrino signal would be feasible also for optically-dark 
supernovae in the Galactic center or in molecular clouds, 
which are heavily obscured by dust. Combined anal- 
ysis of the detected signal along with the thermal (10 
MeV, few seconds long) neutrino burst from the cooling 
neutron star which should occur about 10 hours earlier 
('~ Ri,/vs), could provide important clues about the for- 
mation process of the neutron star and the dynamics of 
the resulting supernova shock. 

For the neutrino parameters inferred from recent at- 
mospheric experiments |20| ] , we expect flavor oscillations 
to produce equal fluxes of muon and tau neutrinos upon 
their arrival to Earth, thus allowing for a "r appear- 
ance" experiment. Since the neutrino signal should coin- 
cide with the 7-ray signal, checking the simultaneity of 
photon and neutrino arrival times will allow to test for 
deviations from Lorentz invariance and from the weak 
equivalence principle. Although the neutrino and pho- 
ton signals are expected to spread over ~ 1 hr, detection 
of a large number of neutrinos will allow to test simul- 
taneity to an accuracy much better then 1 hr, thus pro- 
viding better limits than those derived from supernova 
1987A, where simultaneity was tested to an accuracy of 
only several hours (see [ pl| for a review). 
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